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ABSTRACT 



The invention concerns a mediod foi focusing acoustic 
waves useful for obtaining an image of a field to be observed 
m a dissipative heterogeneous medium (2, 3) around which 
acoustic transducers (Tl-Tn, Pl-I'm) forming an imaging 
network and a target net\\^ork lite method consists in 
following a training step during which pulse responses from 
the medhim are measured between each tiansducei (Ti) of 
the imaging net\\'ork (5) and seveial transducers (Tj) of the 
target network (6); deducing therefrom reference signals to 
be emitted by the transducers of the imaging network to 
produce a focused acoustic pulse in each transducer of the 
tai^et network, then cumulatively, in determining reference 
signals to be emitted to focus an acoustic pulse on prede- 
termioed points in the medium Said reference signals are 
stored and used subsequently to generate an acoustic image 
of the medium . 

23 Claims, 1 Dr awing Sheet 
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METHOD AND N0N-IN\:4SIVE DEVICE FOR transducers are not in dii-ect contact with the abenating 
FOCUSING ACOUSTIC WAVES layer), these impulse responses heing stored in digital form 

viith a certain time sampling which deteimines a number p 
f lELD Of THE DISCLOSURE of frequency components of the impulse response^ with 

5 respective frequencies ok, i being an mdex between 1 and 
This application is a 371 of PCT/FR01/03208 filed on n which designates a tiansducerofthe imaging aiTay^r being 
Oct 1 7, 2001 an index between 1 and m which designates a focusing point 

The present invention relates to noninvasive methods and corresponding to a tiansducei of the target array and k being 
devicesfor focusing acoustic waves, in particular ultrasound an mdex between 1 and p which designates a frequency 
waves. 10 componentj 

(bl2) on the basis of these impulse responses, for each 
BACKGROUND OF THE DISCLOSURE focusing pomt i corresponding to a transducer of the target 

anay, a set of n reference time signals e'i(tjr) is calculated, 
More particular ly, the invention relates to a noninvasive i varying between 1 and n, such that if the aberrating wall 
method for- focusing acoustic waves in a dissipative heteio- 15 were removed in the vicinity of the focusing point r, the 
geneous medium comprising a substantially homogeneous emission of these reference signals by the various ttansduc- 
medium (for example the bmin) surrounded at least partially ers i of the imaging air ay would generate a predetermined 
by a dissipative aberrating layer (for example the skill) signal (for example an acoustic pulse) focused on the 
which generates aberrations in tlie propagation of the acous- focusmg point r, 

tic waves, the acoustic waves being emitted from outside the 20 (b2) a substep of focusing at a number R of predetermined 
aberratmg layer and focused in the substantially homoge- focusing points lying in the substantially homogeneous 
neons medium . medium, with indices q between m+ 1 and m+R, this substep 

The methods of this type which are commonly used do not consisting in determimng for each of these focusing points 
make it possible to obtain good focusing of the acoustic moving step-by-step away from the focusing points 1 to 
waves inside the medium, and, when these methods are used 25 m conespondmg to the transducer s of the target an ay 
in imaging applications, they therefore do not make it reference signals e*i(t,q) to be emitted by the various trans- 
possible to obtain a good resolution and a good image ducers i of the imaging array in or der to generate a pulse 
contrast when the propagation aberrations are significant, for focused on said focusmg point q, the reference signals 
example when echography of the bram is being canied out e'i(t,q) being determined for each focusing pomt q by 
from outside the skull 30 proceeding as follows: 

It is, in particular, an object of the present invention to (|j2i) a fct estimate of e'i{t,q), for i rangmg from 1 to n, 
overcome this drawback jg calculated on the basis of at least one reference signal 

oTTHjmf r^c TTjr: TMcr-T noTmr? ^^(^t'^^)f hemg the iudex of at least one focusing point 

SUMMARY OF THE. DISCLOSURE ^^^^^ focusmg point q for which the reference signal 

, . , ^5 has already been determined, this calculation being per- 

To this end, accordmg to the mvention, a focusmg method ^^^^^ |^ ■ ^ ^ ^^^^^ ^^^^^^-^ ^^^^ ^ 

of the type m question is characterized m that it includes the substantially homogeneous medium (2), 

followmg steps: (b22) the transducers of the imaging ariay are made to 

fa) an initial positionme step diums; which a number t \^ ^ . ^ ^ * • 1 1.* • j jr*t- 

_/ .f ^ r,-,- emit by iteianons, the estmiates previously obtamed of the 

greater than 2 01 acoustic transducers are faxed m predeter- 40 ^ ' . , iv+ \ +t, * 1. . /+ i\ i^^^v .^^.^^..^a 

. .. a. -J f ^.1 1 5 J reierence signals eift.q), then simals s.ft.q) back-scatt^^^^ 

mined positions outside the abenatmg layer, these transduc- . . • 1 ; a- * i.^^„^„,;+t. 

... . 'I 1 . by the dissipative heterogeneous medium are pi eked up with 

ers being m (direct or indirect) contact with said abenatmg transducers, then diese reference signals eXt,q) are 

layer and forming at leas^: ^ , modified for the next iteration in the following way: 
an miagmg an ay which combmes a number n betvi'een 1 

and t of said transducers, 45 (t)-*ai(q) e,- (t-t/q)) 
and a target array which combines a number m between 

1 and t of said transducers (these two arrays may be where the values af(q) and t/q) are a corrective amplitude 

entirely separate, or include certam common transduc- factor and a corrective delay, which are calculated so as to 

ers, or aheinatively each inchide all of the aforemen- maximize a coherence criterion C between said back-scat- 

tioned transducers), 50 signals, said iterations being stopped when the crite- 

(b) a learning step itself compiismg the followmg sub- "o^i ^ reaches a predetermined threshold, 

steps: (b3) the reference signals e'i(t,q) aie stored, at least for q 

0)1) a substep of learning to focus the miagmg array on between m+1 and m+R, 

the target anay, during which substep; (c) and a focusing step during which, for at least one of 

(bll) impulse responses hri(t) of the dissipative hetero- 55 said focusing points q, the transducers of the imagmg array 

geneous medium are determined, respectively between each are made to emit said reference signals e'i(t,q), i bemg an 

transducer iofthe imaging array and a plurahty of focusmg index between 1 and n designating a tr ansducer of the 

points I lymg on the abenating layer in respective corne- imaging anay. 

spondence with transducers of the target anay (this deter- By virtue of these provisions, die propagation aberrations 

mination may be carried out by direct measurement if the 60 of the acoustic waves in the dissipative heter ogeneous 

(lansduceis of the target array are made to emit acoustic medium are overcome and ver>' precise focusing is obtained, 

pulses, or optionally by measurement and calculation if the which may in particular make it possible to obtain reliable 

transducers of the target anay are made to emit acoustic and precise echography ofa field to be observed through the 

signals other than pulses; the values measured and/or cal- abenating layer by hack-scattering, when acoustic waves 

dilated in this way may then optionally be corrected by 65 focused on different points of the field to be observed are 

digital backpropagation in order to simulate transducers successively emitted and the back-scattered acoustic waves 

lying directly in contact with the aberxatmg layer if the are picked up 



us 1,V 

3 

Ibis precise focusing may also be used in applications 
ottier than echography, in paiticuiar: 
Dopplei color imaging, 

elastogiBphic imaging methods, such as the one desctibed 
in document WO'A-00/55 616, 

nonlinear imaging methods ("liaimoiiic imaging"), 
methods of treatment by localized destrxiction of apart oi 

the dissipative heterogeneous medium, in particulai by 

hyperthermia, 

methods foi measuiing optical absorption paiameteis of 
tissues with activation by ultiasound, etc 

hi prefeired embodiments of the invention, one and/or 
other of the following provisions may optionally be hnple- 
mented: 

duiing substep (b 1 1 ), when at least certam transducers (of l s 
the target aiiay and/or the Imaging array) are in contact 
with an intennediate homogeneous medium (foi 
example a gel) which is itself in contact with the 
abeirating iayei, the impulse responses hii(t) are coi- 
rected by digital backpropagation m order to simulate 20 
transducers lying directly in contact with the abenating 
!ayei; 

substep (bl2) itself iacludes the following substeps: 

(bl21) p transfer matrices H(cok)=[Hri(ct)k)] are deter- 
mined, i ranging from 1 to n and i ranging from 1 to m, 25 
where Hri(cok) is the value, at the frequency cok, of the 
Fourier transform of the impulse response hii(t)j 

(bl22) foi each focusing point r corresponding to a 
transducer of the target array, n components Ei(cok,r) are 
determined, i varying between 1 and n, such that F(cok,r)= 30 
H(a)k) E(ct)k,r), where E(a]kj)=[Ei((jL)k^)] is a vector with n 
components, F(cok,i) is a vector with m components Fl(ceik, 
i), 1 vaiyhtg between i and m, these m components Fl(cok,r) 
corresponding to a desired focusing of the acoustic waves at 
the frequency wk on the focusing point r corresponding to a 35 
tr ansducer of the target an ay, 

(bl23) for each focusing point r corresponding to a 
transducer of the target anay, a vector of n time signals 
e(t,r)={ei(t,r)] is deduced therefrom, i varymg between 1 and 
n, where 40 



)1,337 B2 

4 

during substep (bl24), the aberrating wall in the vicinity 
of each focusing point r corresponding to a transducer 
of the target array is assimilated to a filter, which has a 
fimte impulse response and is defined at each frequency 
5 (Dk by an amplitude Gr(ci)k) and a phase substep 
(bl24) itself including the following substeps: 
(bl241) for each fr equency cok, the amplitude Gr((jok) and 
the phase (i)^(ct)k) are calculated on the basis either of the 
signals ei(t,r) or of the vector s E((Dk,r), 
10 (bl242) p corrected transfer matrices H'(cok)=[H'ji(a)k)] 
are calculated, where 



(bl243) for each transducer r of the target array, n 
components E'i(cok,r) are determined, i varying between 1 
and n, such that F(ojk,r)^f^'(Mk) E*(o^k,r), where E'(a)k,i)= 
[Ei(cakjr)] is a vector with n components, F(a3k,r) is a vector 
with m components Vl((£k,t), 1 varying betw^een 1 and m, 
these m components FI(cDkji) correspondkig to a desued 
focusmg of the acoustic waves at the frequency ak on the 
focusing point r conesponding to a transducer of the target 
array, 

(bl244) for each focusing point r corresponding to a 
transducer of the target array, a vector of a reference time 
signals e'(t,i)=[e'i(t,r)] is deduced therefi-om. i varying 
between 1 and n, where 



in complex notation; 
during substep (b2141), the amplitude Gr{o)k) and the 
phase WftJk) are calculated as follows: 



in complex notation, these signals ei(t,r) being adapted so 
that the emission of them respectively by the various trans- 
ducers i of the imaging anay generates an acoustic pulse 
focused on the focusing point r of the target anay, 

(bl24) a substep of correcting the aberrations generated 
by the aberrating layer between the substantially homoge- 
neous medium and each target tr ansducer r, these aberrations 
being estimated on the basis of the measurements carried out 
previously, the aberrations estimated m this way being used 
to calculate said reference time signals e?i(t,r); 
p matrices H"^((Dk) are calculated during substep (bl22), 
respectively by regularization and inversion of the 
transfer matrices H(wk), and the vector E(ajk,r) is 
calculated for each transducer r of the target anay by 
the formula: 

during step (hi 22), the components Fl(a)k,r) of the vector- 
F((i)k,i) coirespondmg to the spatial distribution of the 65 
desired field at the frequency cok are equal to 0 for yi 

and to 1 for H; 



GfiioO = ^ 



where: 

Ei'^ is the complex conjugate value of Ei, 
and AT(rO,r,i)=(d(rO,i)-d(r,i))/c, d(r,i) bemg the dist^ce 
between the transducer i and the focusing point r, and 
d(i€ji) being the distance between the transducer i and a 
particular focusing point rO; 
substep (b]2) itself includes the followmg substeps: 
(bl21) p transfer matrices H(cok)=^[Hri((;ak)] are deter- 
mined, i ranging from 1 to n and r ranging from 1 to m, 
where Flri(cok) is the value, at the frequency tok, of the 
Fourier transform of the impulse response hri(t), 

(bI22') the transfer matrices Fi(a}k) are corrected in order 
to overcome the aberrations generated by the aberrating wall 
m the vicmity of each focusing point r, this correction being 
carried out on the basis of the impulse responses hii(t) 
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determined previously, and coirected transfer matiices 
H'(a)k) are obtained in this way, 

(bl23') foi each focusing point r corresponding to a 
transducer of the target array, n components E'i(wl^r) aie 
determined, i varying between 1 and n, such that F(tokj)= 5 
H'([Dk),B'(ook,r), where E'(tok,i)-[E'i((Dk,i)] is a vector with 
n coioponects, F(cDk,r) is a vector with m components 
Fl(cok,r), 1 varying between 1 and m, these m components 
Fl(a)k,i) corresponding to a desired focusuig of the acoustic 
waves at the frequency wk on the focusing point r corre- 10 
spending to a transducer of the target array, 

(bl24') for each focusing point i conesponding to a 
transducer of the target array, a vector of n tune signals 
e'(t.i)=[e'i(t,r)] is deduced therefrom, i varying between 1 
and n, where 



20 



in complex notation, the signals e'i(t,r) being said reference 

signals; 

p matrices H'"^(cok) are calculated during substep (b 123'), 25 
respectively by regularization and uiversion of the 
transfer matrices ir(a)k), and the vector E^(c}k,r) is 
calculated for each transducer r of the target an ay by 
the foimitla: 

during step (bl23'), the components Fl(cok,r) of the vector 
(cok^r) corresponding to the spatial distribution of the 
desired field at the frequency oik are equal to 0 for yt 
and to 1 for l=i; 35 
during substep (hllT), the aberrating wall in the vicinity 
of each focusing point r corresponding to a transducer 
of the target array is assimilated to a filter ^ which has a 
finite unpulse response and is defined at each JQrequency 
cok by an ainplimde Gr(tok) and a phase (tiX<^k), substep ^ 
(bl22') itself including the following substeps: 
(bl22'l) for each frequency cok, the amplitude Gr(a)k) and 
the phase ^'X^^^) ^® calculated on the basis of the impulse 
responses determined previously, 

(Id 122*2) p corrected transfer matiices H'(ci)k)=[If ji(a)k)] 
are calculated, where 



H'^ri designates the complex conjugate value of Hri, 
and AT(rOjr;i)=(d(rO,i)-d(r,i))/c, d(r,i) bemg the distance 
between the transducer i and the focusing point r, and 
d(rO;i) being the distance between the transducer i and a 
particular focusing point rO; 

durmg step (c), substep (cl) is followed by the following 
substeps: 

(c2) said transducers of the hnaging array are made to 
pick up signals s^(t) back-scattered by the dissipative het- 
erogeneous medium, 

(c3) the reference signal emitted by each transducer of the 
hnaging array is convoluted with the back-scattered signal 
picked up by this transducer, 

(c4) then the convolution products obtabed in this way 
are summed, 

step (c) being repeated for a plurality of points lying in the 

substantially homogeneous medium; 
during substep (b21), the first estimate of each reference 
signal is e'i(t,q)=e'i(ts+9i(q),q0) for each focusing point 
q, qO being the index of a focusing point close to the 
focusmg point q for which the reference signal has 
already been determined, 9i(q) being a delay equal to a 
value 5i(q)/c, where c is the average speed of the 
acoustic waves in the medium, and 6i(q) is equal to a 
difierence between^ on the one hand, a distance 
between the tr ansducer i of the imaging array and the 
focusmg pomt qO, and, on the other hand, a distance 
between the transducer i of the imagmg anay and the 
focusing point q; 
during substep (b2), when at least certain transducer s with 
mdex v of the imaging array are not directly m contact 
with the abenating layer, the corresponding signals 
Q\(t,q) are con*ected by digital backpropagation in 
order to simulate transducers placed ui dhect contact 
with the abertatmg layer; 
duimg substep (b22), the values a^(q) andx^q) are looked 
for to maxkoize the following coherence ciitetion C: 



where ; 
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E 




ir-r, q) 
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n 

n 1: < 


(^■^ S 
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during substep (bl22'i), the amplitude Gr(cok) and the 
phase ^^isy^ are calculated for each fi-equency cok in 
the followmg way: ^5 



M(^k)^-]]{^3LfiMe^^'^''''^^^^ where: 



60 



65 



g;(t,q)=Sj.(tXx) eXt,q), x representing the convolution opera- 
tion, 

and <> represents a time average; 

during substep (b22), the values Tj(q) are calculated by 
maximizing a cross-corr elation function, for transduc- 
ers close to the imaging anay, of the signals g,(t,q) and 

dming substep {b22), the values o^Cq) are calculated so as 
to equahze the maximum amplitude of the fimctions 
g,(t,q) on the index i; 

during substep (b22), the values a^q) and T^q) are 
calculated by carrymg out a cross-corielation, for trans- 
ducers close to (he hnaging array, of the signals gjf(t,q) 

during substep (b22), the values a/q) and t^(q) are 

calculated so as to equalize the maximum amplitude of 

the functions g^-ftq) on the index i; 
substep (b22) relathig to each focusing point q is carried 

out immediately after substep (b2 1) relating to the same 

focusing pomt q; 
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the dissipative heterogeneous medium consiste of the 

biain suirounded by the skull; 
the imaging aiiay and the target aiiay are two separate 
aiiays arranged on either side of the dissipative het- 
erogeneous medium; 
all the ttansdiKers belong both to the imaging aitay and 

to the target aiiay; 
the acoustic waves are ultrasound waves. 
ITie invention fmtheimore relates to a device designed for- 
implementing the method defined above. 

Other characteristics and advantages of the invention will 
become apparent during the following description of one of 
its embodiments^ which is given by way of nonlimiting 
example with refenence to the appended drawing 

BRIEF DESCRIFIION OF IHE DRAWINGS 

In the drawing, the single FIGURE represents an ultra- 
sound imaging ctevice according to one embodiment of the 
invention 

DETAILED DESCRIPIION OF IHE 
DISCLOSURE 

The ultr asound imaging device 1 represented in the draw- 
mg is designed to produce an ultrasound echographic unage 
of a patient* s bmin 2 (at frequencies of, for example^ the 
order of from 1 to 3 MHz) from outside the skull 3, the brain 
2 constituting a substantially homogeneous medium for 
propagation of the acoustic waves and the skull 3 constitut- 
ing a dissipative aberrating layer, so that the overall cranium 
2, 3 constitutes a dissipative heterogeneous medium 

As a variant, the invention could be applicable, in par- 
ticular: 

to acoustic imaging of any other nonhomogeneous dissi- 
pative heterogeneous medium comprising a substan- 
tially homogeneous medium surrounded by a relatively 
thuL dissipative layer generating aberrations in the 
propagation of the ulttasouad waves, 
or to any other method involving at least one fbcusmg 

upon emission into such a medium 
In the example represented in the dr-awing, the hnaging 
device 1 includes a microcomputer 4, or any other device for 
control and^or visualization of ultrasound images, this 
microcomputer conventionally including a keyboard 4a, 
optionally combined witii other control interfaces, and a 
screen 4b making it possible to visualize the images of the 
brain 2 

Ihe imaging device 1 fiu:thetmore includes two arrays 5, 
6 of ultrasound transducers Tl, 12 , In and T% V2 
I'm forming, for example, two linear- banks of transducers 
which are arranged on either side of the user 's skull 3, in 
predetermiaed geometrical positions with respect to one 
another, each transducer bank 5, 6 being brought into contact 55 
with the skull 3 via a layer 7 of gel or the like 

The various transducers Tl, 12 In and I'l. r2 . Vm 
may be controlled dnectly by the microcomputer 4, or 
preferably by a central electronic unit CPU which is con- 
tained, for example, in an electronics rack 8 and is itself 60 
controlled by the microcomputer 4 

Advantageously, each of the transducers II, 12 .. In, 
II, V2, Tm is connected to a sampler, respectively El, 
E2 . En, El, E'2, E'm, and each sampler is itself connected 
to a memory^ respectively Ml, M2 . Mm, Ml, M^2 . . 65 
Mm and to a central unit CI, C2 . Cm, C'2, C2 Cm 
These memories and these central units are in tum con- 
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nected, dnectly or indirectly, to the aforementioned central 
unit CPU, which is iurthermoie connected to at least one 
central memory M. 
Ihe device which has just been described operates as 

^ follows. 

Initially, the two arrays of transducers 5, 6 are fixed on 
either side of the patient's skull 3, in said predetermined 
positions. lb this end, the armys of transducer s 5, 6j respec- 
10 tively referred to as the imaging array and the target anay, 
may be carried by a rigid support such as a hat (not shown) 
arranged around the patient*s head 

The device then follows a learning step lasting a few 
minutes (advantageously from 1 to 3 min), making it pos- 
sible to take account of all the propagation aberrations due 
to the nonhomogeneous nature of the dissipative medium 
formed by the skull 3 and the brain 2 

Dining this learning step, fir stly each of the transducers 
Ilj 12 . li, . Tn, of the imaging aiiay 5 is made to 
successively emit an acoustic pulse, and, for each pulse 
emitted by one of the transducers li of the imaging anay, the 
signal picked up by the transducer's Tl, II T'r, , I'm 
of the target anay 6 is recorded, that is to say the impulse 
response hri(t) of the dissipative heterogeneous medium 
between the transducer i in question of the unaguig array 5 
and each transducer j of the target array 6 

Each impulse response hii(t) is recorded in digital form 
with a certain time sampling which determines a certain 
number p of monochromatic firequency components of the 
impulse response, each corresponding to a frequency mk, k 
being an index between 1 and p 

In the case envisaged here, where at least certain trans- 
ducers of the target anay and/or of the imaging array are not 
directiy in contact with said aberrating layer 3, the impidse 
responses are corrected in order to simulate vktual trans- 
ducer's arranged in contact with said aberrating layer The 
position of the layer- with respect to the transducers may 
optionally be obtained by conventional imaging (ultrasound 
echography, x-ray scanner, MRI, etc) The corrected 
impulse responses are calculated by a known digital back- 
propagation algorithm, described in particular in the follow- 
ing articles: 

**Ultrasonic beam steering through inhomogeneous layers 
with a time reversal mirror'', C, DORME, M. FINK, 
IEEE Transaaiom Ultrasonics, Ferroelectric and Fre- 
quency Control, 43 (1), January 1996, pp 167-175, 
"Focusing and steering through absorbing and aberrating 
layers; Application to ultrasonic propagation through 
the skull" Joiunal of Acoustical Society of America, 
103 (5), May 1998, pp 2403-2410, 
and ''Propagation and backpropagation for ultrasonic 
wavefront design" Liu, D.-L , and Waag, R C IEEE 
Trans, on Ultras Ferro. and Freq Contr. 44(1): 1-1 3 
(1997) 

In what follows. hri(t) will therefore denote the impulse 
responses for (real or virtual) elements lying against the 
aberrating layer; The virtual or real elements lying against 
the aberrating layer 3 will fiirthemiore be referred to below 
as "focusing points" of index r between 1 and m 

When the transducers of the imagmg array 5 emit acoustic 
signals e.(t), these signals generate acoustic signals fi(t) 
expressed as follows at die transducers i of the target anay 
6; 
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Ihe gain factor Gj(cDk) and the phase factoi <t)/Mk) are then 
calculated for each focusing point i: 



where x represents the tiiae convolution opemtion 0 £)C&Ji 

After Fourier transfoina, this equation becomes; V 

^i>A) = ^£(arg(£i(tr>i K)))-aig(E,(^^t. r^-^^^'*^-^^)) 



After Fourier transforna, this equation becomes; 

Fii^)kyH(^ak) £(ojfr}, where: 10 



H(ti)k) is the transfer matrix, of size m*n, between the 
tr ansducer s T'i of the imaging an ay and the transducers 
Tr of the tar^get array: the components Hri(tok) of this 

matrix are the components of the Fourier transforms of where Ei* is the complex conjugate value of Ei. 

the impulse responses hii(t) at the frequency tok, The pair^ {Gj(a)k), ^/o)k)} correspond to the relative 

E(wk) is a vector whose components E,(«k) are the attenuation factor and relative phase shift mtroduced at each 

components of the Fourier tiansfotm of the afbimien- frequency by the portion of the aheriating layer 3 lying 

tioned signals e^t) at the frequency cok, ^S^^nst the focusing point r. They therefore finally charac- 

and F(a)k) is a vector whose components F .(a)k) are the ^^"^^ aberrations introduced by the abeimtmg layer 

components of the Fourier transform of the aforemen- ^onion lying against the target array 

tioned signals f/t) at the frequency wk The aberrations introduced by the abetratmg layer 3 lying 

By iaveititig each transfer matrix H(a)k), h is therefore ^S^^t the target transducers are then diminated in all the 

possibletodeterminethevectorE(cokJ)whichissuitablefor P matrices H(a)kHHji(a)k)] defined above, 

generating, at the focusing pomt r conesponding to tiie 25 ^^'^ ^ ""^^ of transfer- matrices H(o3k)=[Hji 

transducer Tr of the target array, a vector F(wko) all of C^^)] calculated characterizing the propagation between 

whose components ar^ as close as possible to l3ie objective thermagmganay and the target anay ma virtual medium for 

initially fixed (preferably all equal to zero, except for the "^^^^ aberrations lying against the imaging anay 

component with index j conespondmg to the transducer T^j, rcniam: 
which is equal to 1 when the intention is to emit an acoustic 30 
pulse at the focusing point r), by virtue of the relationship: 

where H"'(cok) is the invetse matiix of H((iok) 

H-'(cok) may be calcnlated by singular value decompo- Foi each transducei i oi the tatget airay, n components 

sition, foi example, which makes it possible to regularize the E'i(a)k^) are then determined, i varying between 1 and n, 

inversion of (he matiix H(«)k) such that F((ok,i)=H'(cok) P(ci)k,r), where E'((ok,r)=[E'i(tok, 

Next, by inverse Fourier transform of the various com- « ^^'='«' "'jj'i » components F(ft)k,r) is a vector with 

ponents Ei(K>kj) of the vector E(a)k,j), the various reference ™ components Fl(a)k,r), 1 varying between 1 and m, these Hi 

signals ei(tj) are determined which are suitable for focusing components Fl(a)k,r) correspondmg to a d«ued focusmg o 

an acoustic pulse (or optionally another acoustic signal) at acoustic waves at the fiequency cokonthe focusing pomt 

the focusing point i, when they are emitted by the various ' corresponding to a transducer of the target array 

transducers li of the imaging anay 5 Focusing of the Foreach focusmg pointrcoiresponding to a transduoer of 

imaging anay 5 on each transducer of the target anay 5 is <^ ^'^y' ^ vector of n reference hme signals <^(t^)= 

therefore canied out by an inverse spatiotemporai filter Te'iftr)] is deduced therefrom, t varying between 1 and n, 

Ihe central unit CPU then follows a process of learning 
the aberrations at the target array due to the wall of the skull 

During this process, these aberrations are considered as a r)-J^ Ei{ojk ^) ' 

filter with a finite impulse response i=i 

In the Fourier- domain, this fiher is defined at each 
frequency cok by an amplitude Gr(Q)k) and a phase (|)^(o>k) 

Li order to calculate these coeJB&cients, the phase and the m complex notation 

amplitude ofail the vectors Er are compared. To this end, the 55 These reference signals e'i(tj) are adapted so that the 

first stage is to eliminate the phase shifts introduced by tire emission of them respectively by the various transducers i of 

path differences between the imaging transducers Ti and the the imaging array generates an acoustic pulse focused on the 

vaiious focusing points indexed r Ihis is equivalent to transducer j of the target array in tiie absence of the 

selecting a particulai focusing point rO and introducing a aberrating layer lying against the target anay, 

linear phase shift for the others with the angular frequency: 60 It will be noted that, as a variant, the reference signals 

exp(-jAT(ro,r,l)a)) with AT(rD,rJ)^(d(rO,i)-d(r,i))/c where could be determined in the following way, after having 

d(r,i) is the distance between the transducer i and the determined the impulse responses hri(t) and the p transfer' 

focusing point r, and c is the average speed of the acoustic matrices H(a}k): 

waves m the medium to be imaged, in this case the brain 2. the transfer matrices H(coJ£) are corrected in order to 

Once this correction has been carried out, the differences 55 overcome the aberrations generated by the aberrating 

m amplitude and phase between the vectors Er are attributed wall 3 in the vicinity of each focusing point r, this 

to the abenatmg layer 3 which lie against the target array, correction being carried out on the basis of the impulse 
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responses hriCt) determined previously, and corrected 
transfer matrices E\fX)k) are obtained in this way, 

by inverting tbe matrices H'((i)k), for each focusing pomt 
r corresponding to a transducer of the target array, n 
components E'i(ok,i) are determined, i varying 5 
between 1 and n, such that F(cokii)==H'(cok) E'(Q)k,i), 
where E'{cok^)=^[E'i(cok^)] is a vector with n compo- 
nents, P(cDk.r) is a vector with m components Fl(c]ok,r)j 
1 varying between 1 and m, these m components 
f l((i)kjr) corresponding to a desiied focusing of the to 
acoustic waves at the tirequency cok on the focusing 
point I' cojresponding to a transducer of the target array, 

and for each focusing point i corresponding to a trans- 
ducer' of the target array, a vector of n time signals 
e'(f,r)=[e^i(t,r)] is deduced therefrom^ i varying between ^ ^ 
1 and n, where 



20 



ui complex notation, the signals e'i(t,i) being said reference 
signals 25 

Advantageously, during the calculation oi the matrices 
H'(wk), the aberrating wall in the vicinity of each focusing 
point 1 corresponding to a transducer ot the target array is 
assimilated to a filter, which has a finite impulse response 
and is defined at each frequency cok by an amplitude Gr(e)k) 
and a phase tj)/cDk), which are calculated as follows: 



2 



35 



H*ri designates the complex conjugate value of Hri, 
and AT(rO,r3i)={d(rO,i)-d(r,i))/c, d(i,\) being the distance 
between the transducer i and the focusing point i' and 
d{r03i) being the distance between the transducer i and a 
particular focusing point rO 

p corrected transfer matrices H'(a)k)=[H*ji(cak)] ate then 
calculated, where 



which are used to determine the vectors E'i(Q)k) as explained 
above, and therefore the various reference signals e'(t,r), r 
ranging from I to m 

The centra] unit CPU then learns to focus at a number R 
of predetermined focusing points lymg in the brain 2, with 
mdices q between m+1 and m+R^ this substep consistmg in 
determining for each of these focusing points q, moving 
step-by-step away fix^m the transducers of the target anay, 
reference signals e'i(t,q) to be emitted by the various trans- 
ducers of the imaging anay in order to generate a pulse 
focused on said focusing point q 



The reference signals e'i(t,q) are initially detemiined, for 
each new focusing point q, in the form e'i(t,q)=e'i(t+Gi(q)j 
qO) for each focusing point qO being the mdex of a 
focusing point close to the focusmg point q for which the 
reference sigaal has already been determined, the delay 61(q) 
initially being equal to a value 6i(q)/c, where c is the average 
speed of the acoustic waves in the medium, and 5i(q) is 
equal to a difference between, on the one hand, a distance 
between the transducer i of the imagmg airny and the 
focusmg point qO, and, on the other hand, a distance between 
the transducer i of the unaging array and the focusing point 

q- 

In the event that certam transducers with mdex v of the 
imagmg anay do not lie against the abenating layer, it is 
fiirthermore desirable to correct the reference signals the 
corresponding signals e\{t,q) by digital backpropagation 
from the vhtua] transducers (lymg against the abenating 
layer 3) to the real transducers (separated from said layer 3 
by some gel 7 or the like), in a marmer which is known per 
se, by the method which is the rever se of that described 
above in relation to the impulse responses. 

The transducers of the imaging array are then made to 
emit, by iterations, the estimates obtained for the reference 
signals e'i(t,q), then signals Sf(t,q) back-scattered by the 
dissipative heterogeneous medium are picked up with the 
same transducers 

Next, these reference signals e'i(t,q) are modified for the 
next iteration in the following way: 

ei'CtHa,(q)e,'(t-xM 

where the values Tj(q) and a .(q) are a conective delay and 
a corrective amplitude factoid which ate calculated so as to 
maximize a coherence criterion C between said back-scat- 
tered signals 

Advantageously, this coherence criterion C may be the 
followmg: 
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45 

g,.(t,q)=sXt)xe/(t,q) x representing the convohition opera- 
tion, 

and <> represents a time average 

In this optimization process, the values t^q) may be 
calculated at each iteration so as to maximize a cross- 
correlation function, for transducers close to the imaging 
anay, of the aforementioned signals gf(t,q) and g,-+i(t,q), and 
the values o(^(q) may be calculated so as to equahze the 
maximum amplitude of the functions g(t,q) on the index i 
The reference signals e4(t,q), i ranging from 1 to n, are 
hence optimized so that they produce an acoustic signal 
focused precisely on the focusing point q lymg in the brain. 
Ihis optimization process has already been explained in 
50 more detail by Mallart et al. (The Van Citteit-Zeinike 
theorem in pulse echo measurements^ J.. Acoust Soc. Am. 
90(5), November 1991, pp 2716-2727; Adaptive focusing 
in scattering media through sound speed inhomogeneities: 
the Van Cittert Zemike approach and focusing criterion, J. 
e5 Acoust Soc. Am. 96(6), December 1994, pp. 3721-3732) 
When this optimization is completed for a focusing point 
q, for example after 2 or 3 iterations when the criterion C has 
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reached a piedetennined value (ia particular close to 5^), 
operation proceeds to the next focusmg point q+1, etc 
The reference signals ^i{X,q) obtained in this way are 

stored, for example in the memories Ml-IVln 

Once the learning step is completed^ it is in particular 5 
possible to produce echographic images of the brain 2, 
optionally at a fast rate which may be as high as the speed 
of a staiadaid echogiapb, foi example 20 to 30 images per 
second. In oidei to produce each of these images^ the 
following procedure is adopted foi each focusmg point q lo 
belonging to the field to be observed: 
the transducer s li of the imaging array are respectively 

made to emit said reference signals ei(t,q), 
then said transducers of the imaging anay are made to 
pick up signals si(t) back-scattered by the viscoelastic 15 
medium, 

tlie reference signal ei(t,q) emitted by each transducer of 
the imaging anay is convoluted with the back-scattered 
signal si(t) picked up by this transducer, 

then the convolution products obtained in this way are 20 
summed. 

It will be noted that the various aforementioned opera- 
tions canied out during the learning step or the imaging step 
may either be programmed in the central unit CPU, or all or 
some of them may be performed by specialized circuits 25 

Puithermore, it will also be noted that all the transducers 
Ti, T'l coufd be used to produce the echographic images of 
the biain In this case, the imaging array would be the same 
as the target array, and each of these two arrays would 
comprise all the transducers, the operation described above 30 
then being applied mutatis mutandis 

The invention claimed is: 

1, A noninvasive method for focusing acoustic waves in 
a dissipative heterogeneous medium (2, 3) comprising a 
substantially homogeneous medium (2) suir-ounded at least 
partially by a dissipative aberrating layer (3) which gener- 
ates aber rations in the propagation of the acoustic waves, the 
acoustic waves being emitted fi:om outside the aberrating 
layer (3) and focused in the substantially homogeneous 
medium (2), 

characterized in that it includes the following steps: 

(a) an initial positioning step during which a number t 
greater than 2 oi acoustic transducers (Tl-TUj 
T'l-T'm) are fixed hi predetermined positions outside 45 
the aberrating layer (3), these transducers being in 
contact with said aberrating layer and forming at least: 

an imaging anay (Tl-Tn) which combines a number n 

between 1 and t of said transducers, 
and a target array (PI- I'm) which combines a number m 50 

between 1 and t of said transducer s, 

(b) a learning step itself comprising the following sub- 
steps: 

(bl) a substep of leariung to focus the imaging array on 
the target array, during which substep: 55 
(bll) impulse responses hri(t) of the dissipative hetero- 
geneous medium are determined, respectively between 
each transducer i of the imaging array and a plurality of 
focusing pomts r lying on the abenating layer (3) in 
respective correspondence with tr ansducers of the tat- 60 
get array, these impulse responses being stored in 
digital form with a certain time sampling which deter- 
mines a number p of frequency components of the 
impulse response, with respective firequencies cokj i 
being an index between 1 and n which designates a 65 
transducer of the imaging anay, r being an index 
between 1 and m which designates a focusing point 
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corresponding to a transducer of the target array and k 
being an mdex between ] and p which designates a 
frequency component, 
(bl2) on the basis of these impulse responses, for each 
focusing point r corresponding to a transducer of the 
unaging array, a set of n reference time signals e'i(t,t) 
is calculated, i varying between 1 and n, such that if the 
aberiatmg wall were removed in the vicinity of the 
focusing point r, the emission of these reference signals 
by the various transducers i of the imaging array would 
generate a predetennined signal focused on the focus- 
mg point r, 

(b2) a substep of focusing at a number R of predetermined 
fbcusing points lying in the substantially homogeneous 
medium, with indices q between m+l and m+R, this 
substep consisting in determining for each of these 
focusing pomts q, moving step-by -step away from the 
focusing points 1 to m corresponding to the transducer s 
of the target array, reference signals e^t^q) to be 
emitted by the various transducers i of the unaging 
array in order to generate a pulse focused on said 
focusing point q, the reference signals e'i(t,q) being 
determined for each fbcusing point q by proceeding as 
follows: 

(b21) a first estimate of e'iftq), for i rangmg from 1 to n, 
is calculated on the basis of at least one reference signal 
e'i(tjqO), qO being the index of at least one focusing 
point close to the focusing point q for which the 
reference signal has akeady been determined, this 
calculation being performed by using an average speed 
of the acoustic waves in the substantially homogeneous 
medium (2), 

(b22) the transducers of the imaging array are made to 
emit, by iterations, the estimates previously obtahied of 
the reference signals e'i(t,q), then signals s^(t,q) back- 
scattered by the dissipative heterogeneous medium are 
picked up with the same transducer^;, then these refer- 
ence signals e'i(t,q) are modified for the next iteration 
in the following way: 

where the values aXq) and T^q) ar e a corrective ampli- 
tude factor and a corrective delay, which are calculated 
so as to maximize a coherence criterion C between said 
back-scattered signals, said iterations being stopped 
when the criterion C reaches a predetermined thresh- 
old, 

(b3) the reference signals e^i(t,q) are stored, at least for q 

between m+l and m+R, 
(c) and a focusing step during which, for at least one of 
said focusing points q, the transducers of the imaging 
array are made to emit said reference signals e'Xt,q)j i 
being an index between 1 and n designating a trans- 
ducer- of the imaging array 

2 The method as clauned in claim 1, m which during 
substep (bll), when at least certain transducers (Tl-Tm^ 
Tl-T'm) are in contact with an intermediate heterogeneous 
medium which is itself in contact with the abenating layer' 
the impulse responses hii(t) are corrected by digital back- 
propagation in order to simulate transducers lymg directly in 
contact with the aberrating layer;. 

3 The method as claimed in claim 1 or claim 2, in which 
substep (bi2) itself mcludes the following substeps: 

(bl21) p transfer matrices H(ct)k)=[Hri(wk)] are deter- 
mined, i ranging fiom 1 to n and 1 ranging from 1 to m, 
where Hii(Q)k) is the value, at the frequency cok, of the 
Fourier transform of the impulse response hri(t). 
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(bl22) foi each ibcusing point r coiresponding to a 
tr ansducei of the target an ay, n components Ei(Q)kjf) 
are deteimined, i varymg between 1 and n, such that 

F(a)kj)=H(a)k) E(G>k,i), where E(cok,i)=[Ei{tok,r)] is a 
vectoi with n components, F(cak,r) is a vector with m 
components Fl(tok,i), 1 varying between 1 and m, these 
m components Fl(a)k,r) corresponding to a desired 
focusing of the acoustic waves at the ftequency cuk on 
the focusing point r corresponding to a transducer of 
the target anay, 
(bl23) for each focusing point r coiresponding to a 
transducer of the target array, a vectoi of n time signals 
e(t,r)=[ei(t,r)] is deduced therefrom, i varying between 
1 and n, where 



10 



Fl(o3kj) coiresponding to a desu-ed focusmg of the 
acoustic waves at the frequency cok on the focusing 
point I corresponding to a transducer' of the target arr ay, 
(bI244) for each focusmg point i coiresponding to a 
transduce: of the fai^et array, a vector of n reference 
time signals e'(tr)=[e'i(t,i)] is deduced therefi'om, i 
varying between 1 and n, where 



in complex notation 
15 7 The method as claimed in claim 6, in which during 
substep (b 1 24 1), the amplitude Gr(wk) and the phase ^X^^k) 
are cdculated as follows: 



20 
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in complex notation, these signals ei(t,r) being adapted so 
that the emission of them respectively by the vaiioiJS tmns- 
ducers i of the imaging atiay generates an acoustic pulse 
focused on the focusing point r of the target airay, 
(bl24) a substep of correcting the aberrations generated 
by the ab^rating layer between the substantially homo- 
geneous medium and each target transducer r, these 
abenations being estimated on the basis of the mea- 
smements canied out previously, the aberrations esti- 
mated in this way being used to calculate said reference 
tune signals e^i(t,r) 

4 The method as claimed in claim 3, in. which p matrices 
H"^(ojk) are calculated dming substep (bl22), respectively 
by regular ization and inversion of the transfer matrices 
H(cok), and the vector ECoik^r) is calculated for each trans- 
ducer I of the target array by the formula: 

5 The method as claimed in claim 3 or clarni 4, in which 

during step (bl22), the components Fl((jjk,r) of the vector" 
F(wk,r) corresponding to the spatial distribution of the 
desired field at the fiiequency cok are equal to 0 for Is^i and 
to 1 for H 

6 The method as claimed in any one of claims 3 to 5, in 
which during substep Cbl24), the abenating wall in the 
vicinity of each focusing point i coiresponding to a tr ans- 
ducei of the target array is assimilated to a filter, which has 
a finite impulse response and is defined at each frequency cok 
by an amplitude Gr(Q)k) and a phase (|);.(cok), substep (bl24) 
itself including the following substeps: 

(bl241) for each frequency ok, the amplitude Gr(cDk) and 
the phase ij),.(cok) are calculated on the basis either of the 
signals ei(t,r) or of the vectors E(cok,i), 

(bl242) p corrected transfer matrices H'(o]k)=[H'jiCcok)] 
are calculated, where 
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where: 

Ei* is the complex conjugate value of Ei, 

and AT(rO,r5i)=(d(iO,i)-d(r,i)yc, d(r,i) being the distance 
between the transducer i and the focusing point r, and 
d(rO,i) being the distance between the transducer i and 
a particular focusing point rO 

8 The method as claimed in claim 1 or claim 2, in which 
substep (bl2) itself includes the following substeps: 

(bl21) p transfer matrices H(oQk)=[Hri(a)k)] are deter- 
mined, i ranging from 1 to n and r ranging from 1 to m, 
where Hri(cDk) is the value, at the f irequency cok, of the 
Fourier transform of the impulse response hri(t), 

(bl22') the transfer matrices H(tok) are corrected in order 
to overcome the aberrations geneiated by the aberrating 
wall m the vicmity of each focusing point i, this 
correction being canied out on the basis of the impulse 
responses liri(t) determined previously, and corrected 
transfer matrices H'(cok) are obtained in this way, 

(bl23') for each focusmg point i coiresponding to a 
transducer of the target array, n components E'i(cok,r) 
are deterrmned, i vai^'mg between 1 and n, such that 
F(a)k,r)=H'(cok)EXci)k,r), where RCojk,i)=IE'i(o}k,r)] is 
a vector with n components, F(cok,r) is a vectoi with m 
components Fl(ci3k,r), 1 varying between 1 and m, these 
m components F'l(c3kjr) coiresponding to a desired 
focusing of the acoustic waves at the frequency cok on 
the focusmg pouit r corresponding to a transducer of 
the target array, 

(bl24^) for each focusing point r corresponding to a 
transducer of the target array, a vector of n time signals 
e'(t,r)=[e'i(t,r)] is deduced therefrom, i varying between 
1 and n, where 



(bl243) for each transducer i of the target array, n 
components E'i(cok,f) are determined, i varying 
between 1 and n, such that F(tfik,r)=H*(cok).E'(cok,r), 
where E'(cok,r)=[E'i(ft)k,r)] is a vector with n compo- 
nents, F(cj)k5r) is a vector with m components Fl(a)kjr), 
1 varj^mg between 1 and m, these m components 
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in complex notation, the signals e^i(t,i) being said reference 
signals 

9 The method as daimed in claim 8, in which p matrices 
H'"\a)k) are calculated during substep (bl23'), respectively 
by regularization and inversion of the transfer matrices 5 
H'(wk), and the vectoi E^(o)k,r ) is calculated foi each tians- 
ducer r of the target airay by the fonnnla: 

10 The method as claimed in claim 8 oi claim 9, in which 

during step (bl23'), the components Fl(iok,r) of the vectoi 
F(wkji) conespondbg to the spatial distribution of the 
desir ed field at the frequency ok are equal to 0 for Ut and 
to 1 for 1=1 

11 . The method as claimed in any one of claims 8 to 10, ^5 
m which during substep (bl22'), the aberrating wall in the 
vicinity of each focusing point 1 corresponding to a trans- 
ducer oi the target array is assimilated to a filtei' which has 
a finite impulse response and is defined at each frequency cok 
by an amplitude Gr(cok) and a phase ^X^^i substep (bl22*) 20 
itself including the following substeps: 

(bl22'l) for each frequency a)k, the amplitude Gr(cok) and 
the phase ^Xpk) are calculated on the basis of the 
impulse responses determined previously, 
(bl22'2) p corrected transfer matrices H'(mk)=[H'ji(eok)] 25 
are calculated, where 



each reference signal is e'i(t,q)=e'i(ts+6i(q),q0) for each 
focusing point q, qO being the index of a focusing point close 
to the focusing point q for which the reference signal has 
already been determined, 6i(q) being a delay equal to a value 
6i(q)/c, where c is the average speed of the acoustic waves 
in the medium, and Si(q) is equal to a difibrence between, on 
the one hand, a distance between the transducer i of the 
imaging array and the focusing point qO, and, on the other 
hand, a distance between the transducer i of the imaging 
aiiay and the focusing point q 

15 The method as claimed in any one of the preceding 
claims, in which dming substep (b2), when at least certain 
tr ansducer s with index v of the imaging anay are not directly 
in contact with the abenating layerj the corresponding 
signals e'^(t,q) are corrected by digital backpropagation m 
order to simulate nansducers placed in direct contact with 
the aberrating layer: 

16 The method as claimed in any one of the preceding 
claims, in which dming substep (b22), the values a,(q) and 
T^(q) are looked for to maximize the following coherence 
criterion C: 



c = - 
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12 The method as claimed in claim 11, m which during 
substep (bl22'l), the amplitude Gr (a)k) and the phase (jj^wk) 
are calculated foi each frequency ok in the foUowing way: 



40 
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H*ri designates the complex conjugate value of Hri, 
and At(r04,i)=(d(i-0,i)-d(r,i))/Cj d(r,i) being the distance 

between the transducer i and the focusing point r, and 

d(r 0,i) being the distance between the transducer i and 

a particular focusing point rO 
13 The method as claimed in any one of tlie preceding 
claims, in which dming step (c), substep (cl) is followed by 
the following substeps: 55 
(c2) said transducers of the hnaging anay are made to 

pick up signals s/t) back-scattered by the dissipative 

heterogeneous medium, 
(c3) the reference signal emitted by each transducer of the 

imaging array is convoluted with the back-scattered 50 

signal picked up by this transducer, 
(c4) then the convolution products obtained in this way 

are summed, 

step (c) being repeated for a plmality of points lymg in the 
substantially homogeneous medium 65 

14 . The method as claimed in any one of the preceding 
cJaims, in which during substep (b21), the first estimate of 



(=1 



g^(t,q)=s,{t)xe'^(t,q), x representing the convolution opera- 
tion, 

and <> represents a time average. 

17 The method as claimed in claim 16, in which dming 
substep (b22), the values %,{q) are calculated by maximizing 
a cross-coirelation flinction, for transducers close to the 
imaging anay, of the signals g.(t,q) and g,^i(f,q) 

18. The method as claimed in claim 16 or claim 17, in 
which dming substep (b22), the values oL^(q) are calculated 
so as to equalize the maximum amplitude of the fiinctions 
g,.(t,q) on the index i. 

19 The method as claimed in any one of the preceding 
claims, in which substep (b22) relating to each focusing 
point q is carried out innnediately after substep (b21) 
relating to the same focusing point q 

20 . The method as claimed in any one of the preceding 
claims, in which the dissipative heterogeneous medium 
consists of the brain surrounded by the skull 

21 The method as claimed in any one of the parecedkig 
claims, in which: 

either the imaging array and the target anay are two 
separate ariays ananged on either side of the dissipa- 
tive heterogeneous medium, 

or all the transducers belong both fo the imaging array and 
to the target array. 

22 The me&od as claimed in any one of the preceding 
claims, in which the acoustic waves are ultrasound waves 

23 . A device (1) designed for cairymg out a method as 
claimed in any one of the preceding claims, this device 
including a number t greater than 2 of acoustic transducers 
(Tl-Tn, ri-Tm) intended to be fixed in predetermined 
positions outside the aberrating layer (3), these transducers 
being controlled by at least one central electronic unit (CPU) 
and fonning at least: 

an imaging array (Tl-Tn) which combines a number n 
between 1 and t of sard transducers. 
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and a target anay (T'l-rm) which combines a number m 
between 1 and t of said tiansducers, the central elec- 
tronic unit bemg designed to follow the following 

steps: 

(b) a learning step itself compiising the following 5 
substeps: 

(bl) a snbstep of leaining to focus the imaging array on 
the taiget array, during which substep: 

(bll) impulse responses hri(t) of the dissipative het- 
erogeneous medium are determined, respectively lo 
between each iriansducer' i of the imagmg array and 
a plurality of focusing points r lying on the abenat- 
ing layer in respective correspondence with trans- 
ducers of the target array, these impulse responses 
being stored iu digital form with a certain time 15 
sampling wluch determines a number p of fi'equency 
components of the impulse response, with respective 
frequencies ook, i being an index between 1 and n 
which designates a transducei of the knaging array, 
r being an mdex between 1 and m which designates 20 
a focusing point corresponding to a transducer of the 
target aiiay and k being an index between 1 and p 
which designates a frequency component, 

(bl2) on the basis of these impulse responses, for each 
focusing pomt r corresponding to a transducer of the 25 
imaging array, a set of n reference time sigfiak e'i(t^) 
is calculated, i varying between 1 and n, such that 11 
the aberrating wall were removed in the vicinity of 
the focusing point r, the emission of these refer ence 
signals by the various transducers i of the knaging 30 
array would generate an acoustic pulse focused on 
the focusing point i, 

(b2) a substep of focusmg at a number R of predeter- 
mined focusing points lymg m the substantially 
homogeneous medium, with indices q between m+1 35 
and m-^R, this substep consisting in determining for 
each of these focusing points q, moving step-by-step 
away fi:om the focusing points 1 to m corresponding 
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to the transducers of the target array, reference 
signals e'i(t,q) to be emitted by the various transduc- 
ers i of the imaging anay in order' to generate a 
predetermined signal focused on said focusing point 
q, the reference signals e^i(t,q) being determined for 
each focusing point q by proceeding as follows: 

(b21) a first estimate of e4(t,q), for i ranging from 1 to 
4, is calculated on the basis of at least one reference 
signal e'i(tjqO), qO being the index of at least one 
focusing point close to the focusmg point q for which 
the reference signal has already been determined, 
this calculation being performed by using an average 
speed of the acoustic waves hi the substantially 
homogeneous medium (2), 

(b22) the transducers of the unaging array are made to 
emit, by iterations, the estimates pieviously obtained 
of the reference signals e'i(t,q), then signals s^ftq) 
back-scattered by the dissipative heteiogeneous 
medium are picked up with the same transducers, 
then these reference signals e*i(tq) are modified for 
the next iteration in the followmg way; 

s;{t)^iilq).e;{t-xlq)) 

where the values a,(q) ajid i.(q) are a corrective 
amplitude factor and a corrective delay, which are 
calculated so as to maximize a coherence criterion C 
between said back-scattered signals, said iterations 
bemg stopped when the criterion C reaches a prede- 
teimined threshold, 

(b3) the reference signals e'i(t,q) are stored, at least for 
q between m+l and m+R, 

(c) and a focusing step during which, for at least one of 
said focusing points q, the transducers of the imaging 
array are made to emit said refetsnce signals e'i(t,q), 
i being an index between 1 and n designatmg a 
transducer of the imaging array 



